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ABSTRACT 

The existence of correlations between nuclear properties of galaxies, such as the mass of 
their central black holes, and larger scale features, like the bulge mass and luminosity, 
represent a fundamental constraint on galaxy evolution. Although the actual reasons 
for these relations have not yet been identified, it is widely believed that they could 
stem from a connection between the processes that lead to black hole growth and 
stellar mass assembly. The problem of understanding how the processes of nuclear 
activity and star formation can affect each other became known to the literature as the 
Starburst-AGN connection. Despite years of investigation, the physical mechanisms 
which lie at the basis of this relation are known only in part. In this work, we analyze 
the problem of star formation and nuclear activity in a large sample of galaxies. We 
study the relations between the properties of the nuclear environments and of their 
host galaxies. We find that the mass of the stellar component within the galaxies of our 
sample is a critical parameter, that we have to consider in an evolutionary sequence, 
which provides further insight in the connection between AGN and star formation 
processes. 

Key words: galaxies: active - galaxies: evolution - galaxies: ISM - galaxies: nuclei 
- galaxies: stellar content 



1 INTRODUCTION 

Recent statistical investigations pointed out that the nuclear 
regions of approximately 2 per cent of nea rby galaxies show 
emission from ionized gas l|Taniguchi| [2003V The presence of 
this component is mostly explained by processes related ei- 
ther to intense star formation or to non-thermal activity in 
the nucleus. In principle, it is possible to distinguish the sig- 
nature of the two different possibili ties from the spectrum 
emitted by the ionized ga s (see e. g. IVeilleux fc Osterbrockl 
Il987l : iKewlev et all 120061 ). From a practical point of view, 
however, the task is made difficult because of the trend, ob- 
served in large samples of sources, to settle on to a smooth 
sequence, rather than forming distinct distributions. The 
problem is partially due to the difficulties related to the 
spectroscopic measurements, which lie at the basis of the 
technique, but the existence of a population of transition 
objects, characterized by intermediate properties, also plays 
an important role. 
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There are in fact numerous works where the relation- 
ship between nuclear activity and star formation processes 
in the central regions of galaxies is addre ssed, mainly 
takin g i nto account Seyfert 2 galaxies (e. g. | lvanov et al.l 



20001: 



Gonzalez Dclgado, Hcckman & Lcithcrcr 
Gu, Dultzin-Hacyan & Dc Diego 2001; Joguet et al 



2001 



2001; 



20041 ). 



Storchi-Bergmann et alj 200fj : ICid Fernandes et alj 
The properties of Type 1 objects are clearly harder to 
understand, because the light, coming directly from the 
nuclear source, usually outshines the surrounding regions. 
This problem leaves open the question whether the sources 
are intrinsically different or not, which is a concern both 
for black hole-host galaxy relationships and for AGN unifi- 
cation theories. At present, evidence for recen t nuclear star 
form ation exists also for Type 1 AGNs (e. g. iDavies et al.l 
2007) and it is believed that there is little or no difference 
with respect to Type 2 objects under this point of view, 
but the debate is still in progress. 

As a result, though some kind of connection between 
nuclear activity and star formation must exist, its phys- 
ical role is not yet clarified. Several theories were devel- 
oped, mainly invoking effects that could lead from nu- 
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clear star formation to AGNs and then evaluati ng the en- 
erget i c feedback on the surrounding medium (Wccdman 



1983; 



Norman fc Scovilld 11988); IScoyille Sz Norman 
Heckman et al.ll 19891 ; lTaniguchilll999T lEbisuzaki et all 



1988 



2001 



Mouri fc Taniguchj |2002l . However, the possibility that the 



relation may change, depending on the host galaxy proper- 
ties, should also be taken into account. Indeed, it has been 
pointed out that the high velocity outflows, produced by the 
most massive stars of a young stellar population, are very 
unlikely to feed nuclear activity and they can even sweep 
the interstellar medium of low mass galaxies, preventing the 
onset of further star formation or nuclear activity. In very 
massive galaxies, on the other hand, the large amount of 
gas trapped in the gravitational field could serve as fuel for 
a powerful AGN, which is again expected to s top the forma- 
tion of new stars (M onaco fc Fontanodl2005l ; iFerugl io et al.l 
|2010| . and references therein). Such triggers may explain the 
reason why star formation is currently observed over a typ- 
ically limited range of medium mass galaxies. 

In this work, we study the properties of the circum- 
nuclear regions of galaxies, whose spectra indicate the exis- 
tence of an interstellar medium (ISM), which is ionized ei- 
ther by hot young stars or by non-thermal activity of the nu- 
cleus. After selecting our targets, we carry out an extensive 
analysis of their stellar populations, their chemical compo- 
sitions, and their ionization sources. Comparing such prop- 
erties with those of quiescent galaxies, we find additional 
evidence of close relationships between nuclear activity and 
star formation processes, together with several sequences, 
indicating a possible evolutionary connection between the 
different classes. 

Our work is organized as follows: in §2 we describe the 
selection of our sample and its basilar properties; §3 intro- 
duces the data analysis and the strategies that we applied to 
perform our measurements; in §4 we present a discussion of 
our results and of their interpretation, while our conclusions 
are summarized in §5. 



2 SAMPLE SELECTION AND DESCRIPTION 

To investigate the relationship between star forming galax- 
ies (SFG) and AGNs, our work required the selection of 
a wide sample of emission line galaxies. T he Sloan Digital 
Sky S urvey Data Release 7 (SDSS DR7, see lAbazaiian et al.l 
2009) provided the ideal archive of spectroscopic data for 
our purposes. Here we looked for the spectra of galax- 
ies with ionized gas components in the nucleus. The clas- 
sification of these objects according to their ionization 
source led to the introduction of a diagnostic technique, 
which could successfully manage line emitting sources, in- 
cluding both Type 1 and Type 2 AGNs. Our choice fell 
on the diagram plotting log(F [O ii]3727/-qoiii]5007) against 
Iog(F [ oi]630o/-P 1 [olll]5007) (hereafter the O123 diagram). This 
instrument is a combination of two diagnostic ratios that 
express, respectively, the ionization degree and the hard- 
ness of the ionizing spectrum at large radii from the source 
|Kewlev et al ] |200d ; IShields fc Filippenkd[l990h . It provides 
the clear advantages of being compatible with all types of 
AGN spectra, to be nearly independent from the chemical 
composition of gas, and to be free from stellar absorption 
contaminations, which are relevant for most recombination 
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Figure 1. The sample of 119223 emission line galaxies plotted 
on the O123 diagnostic diagram. The red contours show the den- 
sity distribution of sources in levels spaced by 100 objects per 
0.05 dex square bin on the diagram. The blue crosses mark the 
bins corresponding to local minima in the distribution, while the 
blue curve is the borderline defined in Eq. (1), below which 90 
per cent of the objects are AGN powered sources. 



lines. On the other hand, it may suffer from differential ex- 
tinction effects, while the [O I] emission line is often weak 
and difficult to measure in noisy spectra. For this reason we 
collected the spectra of galaxies with detected emission lines 
of [O III] A5007, [O II] AA3727, 3729 and [O I] A6300, requir- 
ing a S/N > 3 for [O I]. We extracted 119226 spectra that 
we plot in the O123 diagram of Fig. 1. 

Comparing the O123 diagnostic diagram with the 
Veilleux-Oste r brock (VO) classification, updated by 
iKewlev et all (|2006h . it can be realized that the observed 
sequences are generally populated by different sources: the 
upper sequence is mostly consisting of star forming galaxies, 
the lower of active nuclei, and the overlapping region of LIN- 
ERs and transition sources. Subdividing the plot in bins of 
0.1 dex in the range -1.2 < log(f] O ii]3727Aqoiii]5007) < 0.3, 
where the separation between the sequences is clearer, 
we computed the number of sources within 0.1 dex wide 
bins of log(F[oi]630o/^[oiii]5007)- We, thus, built histograms 
showing the two-peaked source distributions and we located 
the minimum for each bin. The locations of such minimum 
points were subsequently interpolated by a polynomial 
function, derived by means of a least squares regression, in 
the form of: 



log 



Hon] 3727 
7 [0 III] 5007 



0.20 - 0.25 log 



[O I] 6300 
F[0 III] 5007 



-0.39 



log 



-HOI] 6300 
7 [0 III] 5007 



(1) 



By taking into account only the objects below the curve, 
we isolated a sample of approximately 16000 objects that 
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Figure 2. Examples of stellar continuum spectral synthesis for the case of a Seyfert 2 galaxy (left panel), a normal galaxy (central 
panel) and a star forming galaxy (right panel). The continua were normalized with respect to the flux at 4100 A and the plots show 
the observed spectra in black, the synthetic stellar continua in red and the fitting residuals in cyan. We can easily note the Balmer line 
absorption features that may significantly affect the ratio of important diagnostic lines, especially in the case of S2Gs. 



consists of 60 per cent Seyfert galaxies, including both broad 
and narrow line sources, 31 per cent composite Seyfert - H II 
galaxies, 4 per cent LINERs, and 5 per cent SFGs. 

A detailed classification of the narrow line emitting 
objects, identified by means of a line profile fitting code 
described in a companion paper (GGFIT, IVaona et all 
lacceptedh . can be performed by application of the classic 
VO diagnostic diagrams. Before proceeding, however, we 
have to account for the stellar contributions, which signif- 
icantly affect the intensities of recombination lines like Ha 
and H/3. These contributions were estimated with a spec- 
tral synthesis technique, b ased on the STARLIGHT code 
i|Cid Fernandes et al.ll2005l ). using a set of 92 simple stellar 
population spectra (SSP), obtained combining 23 ages (from 
10 6 yr up to 13 - 10 9 yr) and 4 metallicitie s (Z = 0.004, 0.008, 
0.020 and 0.050) from the library of iBruzual fc Charlotl 
( 2003). In order to fit the observed spectra, we had to correct 
them for the extinction arising in our own Galaxy, accord- 
ing to the extinction map provided by the NED servic^E to 
remove the cosmological redshift, and to re-sample the data 
in a fixed wavelength dispersion. We performed these oper- 
ations by means of the IRAF tasks deredden, newredshift, 
and dispcor. Such procedures were meant to derive accurate 
models of the spectral continua and absorption lines, thus 
producing residuals of pure emission line spectra, provided 
that the S/N of original data was high enough. A threshold 
value of S/N > 10 at 5500 A granted the possibility to apply 
the fitting procedures with reliable results. 

Since we are interested in the properties of the circum- 
nuclear regions of the galaxies, we restrict the sample to 
objects in the redshift range of 0.04 < z < 0.08, where 
the SDSS fiber aperture covers a region of 2.4 — 4.8 kpc in 
radius. In this study we considered only Seyfert 2 galax- 
ies (S2G), collecting 2153 objects. A set of 1302 SFGs was 
identified on the diagnostic diagrams, by applying the con- 
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straints l|Rafanelli et al.ll2009l ): 

log([0 III]/H/3) < 0.61/[log([N II] /Ha)] + 1.3 (2a) 

log([0 III]/H/3) < 0.72/[log([S II]/Ha)] + 1.3 (2b) 

log([0 HI]/H/3) < 0.73/[log([O I]/Ha)] + 1.33 (2c) 

We finally collected an additional set of galaxy spectra, with- 
out detected emission lines, selecting 3000 objects, above the 
adopted S/N > 10 threshold at 5500 A, that were observed 
in the same redshift range. These formed our Normal Galaxy 
sample (NG), which we compared to the line emitting ob- 
jects. More general properties of the AGN cla ss, including 
broad line emitting sources, will be dealt with in lVaona et al.l 
lacceptedh . 

3 DATA ANALYSIS 
3.1 Stellar populations 

The circum-nuclear stellar populations of our galaxy sample 
contribute to the observed spectra with various components. 
Depending on the object distance, different fractions of the 
target galaxy lie in the region covered by the spectrograph 
fibers. As a result, we observe the composite spectrum of in- 
tegrated stellar populations, with their characteristic contin- 
uum and absorption lines, overlapped to the emission lines, 
produced by the diffuse ionized gas. To apply the classifi- 
cation procedures, described above, and to investigate the 
physical properties of the sample, we need to distinguish 
between these components. For this purpose, we reproduced 
the stellar populations, which originate the spectra of all the 
selected galaxies, using STARLIGHT. The models were com- 
puted applying linear combinations of the 92 synthetic spec- 
tra of SSPs, with homogeneous age and chemical composi - 
tion, extracted from the library of lBruzual fc Charlotl (|2003l ) 
and sampled at the same spectral resolution of data. Some 
model examples, for the spectra of NGs, S2Gs and SFGs, 
are given in Fig. 2. After masking the wavelength ranges 
affected by strong emission lines, we modeled the spectra 
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Table 1. Single cloud model parameter space 
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Figure 3. Average normalized spectra of the stellar component 
in the circum- nuclear regions of NGs (black line), S2Gs (red line) 
and SFGs (blue line). 



with STARLIGHT, thus estimating the weight of each SSP 
contribution to the integrated signal and reproducing pure 
stellar spectra. 

Applying this analysis to the whole sample, we observe 
that, on average, the properties of stellar populations differ 
significantly in the nuclei of distinct galaxy classes. Fig. 3 il- 
lustrates this general result in a plot of the average spectra of 
SFGs, S2Gs, and NGs, for comparison, after normalization 
to the specific flux at 5500 A. According to expectations, 
the spectra of SFGs are clearly different, with respect to the 
other classes, due to the presence of young hot stars, which 
dominate the emission, yielding a blue peaked spectral en- 
ergy distribution (SED). The spectra of NGs and S2Gs, on 
the other hand, are remarkably similar and characteristic 
of an evolved stellar population. In spite of the comparable 
SEDs, however, S2Gs show on average a 4000 A break that 
is clearly shallower than the corresponding feature of NGs. 
It becomes therefore evident that an excess of hot stars, 
specifically belonging to the spectral class A, populates the 
circum-nuclear region of active galaxies. 

In the following sections we shall discuss how the av- 
erage properties detected in the sample affect the central 
environment of the objects under investigation. The focus 
of our analysis will be centered on the region covered by the 
2.4 — 4.8 kpc aperture of the SDSS fiber including the galaxy 
nuclei. We dub such region the circum-nuclear part of the 
galaxy, while we talk about nuclear processes to refer to the 
unresolved source of ionizing radiation in the nucleus. 



3.2 Gas components 

Once the stellar contributions in our sample of spectra had 
been reproduced by means of proper templates, we turned 
our attention to the analysis of the diffuse ISM. This task re- 
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quires different approaches for the various classes of galactic 
nuclei we previously distinguished. Indeed, the physical pro- 
cesses, which are relevant to ionized gas in presence of either 
an AGN or young stars, represent far extremes of the ion- 
ization phenomenology and they involve quite different gas 
phases. Each class must, therefore, be dealt with by means 
of the most appropriate technique, though tests performed 
by comparing the results of multiple methods on the same 
objects have also been taken into account. 



3.2.1 Photoionization models 

Modeling the physics of ionized gas emitti n g reg ions in 
AGNs is a great challenge. Following iNetzerl (120081 ) . there 
are five main aspects to be taken into account for a correct 
representation of the ionized gas: photoionization and radia- 
tive recombination, thermal balance, ionizing spectrum, gas 
chemical composition, and cloud or filament distribution. 
In addition, the thermal balance has to consider mechani- 
cal heating and the role of dust, which is always present in 
the narrow line regions (NLR) of AGNs. We have excluded 
from our analysis gas kinematics, winds, and gas confine- 
ment mechanisms. Expl oiting photoionizat ion codes, such 
as CLOUDY v06.02 (see lFerland et al.lll998T ). we are able to 
solve numerically the ionization and thermal structure of a 
single cloud. However, we can also explore different param- 
eter spaces and build more sophisticated models, such as 
composite models with more clouds and various geometries, 
as well as time dependent models. In our calculations we 
assumed open geometry. The input parameters, that we are 
considering, include the ionizing spectrum (a non-thermal 
power law, in our case), the density n e , the ionized column 
density N c (fi II), the ionization parameter U, and the metal 
abundances Z/Zq. The ionization parameter expresses the 
ionization degree of photo-ionized gas and it is defined as 
the ratio of the ionizing photon density and the gas density, 
according to the relation: 



U = 



Q(H) 
47rr 2 criH 



(3) 



In this work we exploited three kinds of photoioniza- 
tion models: a set of single cloud models and two of com- 
posite models, based on a combination of different single 
cloud models. The calculations were developed to explore 
th e parameter s pace, according to the procedure described 
in IVaonal (|201p| ). Actually, in the literature there are two 
fundamental approaches regarding composite models, either 
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cons idering a dense cloud inside a me dium with low den- 
sity l|Binette fc Storchi-Bergmannl 1996ft. o r combining sepa- 
rately two clouds ( Komossa fc Schulzl 1997ft . Since our aim is 
to estimate the mean physical parameters, which character- 
ize the NLRs, using spectra in the visible range, a particular 
attention is put on the metallicity, because this parameter 
is important in the investigation of the nature of gas and of 
its origin. We generated a great number of models, in order 
to fit all the observed lines, subsequently applying a \ 2 test 
to compare our simulations with the observed fluxes. 

In order to explore the 5-dimensional parameter space 
we built a set of single cloud models with the parame- 
ter values reported in Table 1. The ionizing spectra were 
assumed to be power laws with F v cx v a in the range 
10 am - 50keV, a cut-off at low energies {v 2 ' 5 ) and an- 
other at high energies (v~ 2 ). This choice is in agree- 
ment with the results given b y iKomossa fc Schulzl | | 1997ft . 
iGroves. Dopita fc Sutherland! (I 2004ft and supported by X- 
rav observations fsee lMainieri et al.ll2007ft . The basic model 
is conceived to evaluate the ionizing photon flux across a 
slab with a fixed density, column density and dust to gas 
ratio (D/G). The distance of the source does not need to 
be determined, since the ionization parameter is fixed. The 
intrinsic spectrum is calculated using CLOUDY and the cal- 
culation stops when either the fixed required column density 
has been reached, or the temperature in the simulation has 
fallen below 4000 K. We assume that the intrinsic spectrum 
is reddened by the ISM, so that, when the models are com- 
pared to the observed spectra, a correction for extinction 
must be taken into account. 

When the ionizing spectrum is a non-thermal power 
law, the ionization structure is very complex. The partially 
ionized region is so large that different ionization states are 
mixed and the solution of the ionic abundances, which make 
up the total abundances, gets quite complicated. An appro- 
priate choice of the chemical abundance sets is a critical step 
in the creation of the models, because they are strongly con- 
nected with another fundamental aspect of the simulations, 
specifically the role of dust and its composition. Once a ref- 
erence sample of different abundances is assumed (usually 
the solar set of values), we need a law describing the change 
in the chemical composition that will be used for the ex- 
ploration of different metallicities. Observations of the ion- 
ized gas chemical abundances in H II regions and starburst 
galaxies show that all elements, except Nitrogen and Helium, 
in first approximation are found in fixed proportions (e. 
Masegosa, Moles & Campos-Aguilar 1994; Ivan Zee et al.l 



Table 2. Abundances and depletion factors for D/G = 1.00 



1 19981 : Peimbert fc P cimbcrt 2010, and references therein) 



Usually, the H II metallicities are expressed in terms of Oxy- 
gen abundance, because in the visible spectra it is possible 
to measure three different ionization states of Oxygen, which 
is the most abundant element after Hydrogen and Helium. 
Then, we can assume that all elements scale directly with O 
abundance, except for N and He that must be treated sepa- 
rately. The abundance of N, as well as He, must be corrected 
for se condary production (see e. g. IVila Costas fc Edmund! 
1 19931 ): 

X(N) _ X(0) 



X(H) X(H) 
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= 0.0737 + 0.0293^. 
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The composition and quantity of dust are also sources 
of uncertainties. The set of abundances needs to take into 
account both gas and dust composition in order to main- 
tain the total assumed element abundances. As for the 
abundances of the elements, whose reference values are as- 
sumed to be similar to the solar values, the default values 
of abundance and composition for dust are taken from the 
Galactic ISM. Chosen the features of the dust, the abun- 
dances of the gas component are calculated by subtract- 
ing the dust abundances from the total. With the excep- 
tion of the noble gases, all abundances must be corrected 
for depletion into dust. The determination of these correc- 
tions is a main issue of modeling, since the only possible 
assumptions are based on the observational properties of 
the Galactic ISM. Another important assumption involves 
D/G, which we again express in terms of its value relative to 
the local ISM. By changing the relative internal dust con- 
tent as a function of the metallicity, the depletion factors 
must also be changed. A method to estimate the depletion 
fa ctors with the variat ion of the dust content is presented 
m iBinette et all (|l99St ). It has been suggested that the de- 
pletion into gra ins could be c onstant in many directions of 
the Galaxy fsee lVladildl2002 | ). so we adopted the depletio n 
factors reported in lGroves. Heckman fc Kauffmannl (| 2006ft . 
as reference values, and a linear law for scaling the dust 
abundance with metallicity. Using the conventional logarith- 
mic notation, the relation between depletion factor and gas 
abundance is given by: 



de P l(X) = log(X/H) 9as - log(X/H) tot 



(•») 



In this work, we assume log(A/0) = const, for 
all chemical elements, with the exception o f N and He. 
For th ese elements, we follow indications by IGroves et al.l 
( 2004), employing a linear combination of the pri- 
mary and secondary components of N with the require- 
ment to match the adopted sola r abundance patterns 
l|Asplund. Grevesse fc Sauvall l 2005). We focused our atten- 
tion on the high metallicity models, because nuclear gas and 
Seyfert galaxies ver y rarely present sub-solar metallicities 
|Groves et alj r2006). For each adopted D/G we changed 
the depletion factor set in order to retain the total abun- 
dances. The N depletion factor must change with metallic- 
ity because, in agreement with our assumptions, if the dust 
increases with metallicity, there must be an excess of N in 
gas form due to Eq. (4a). In Table 2 we report the deple- 
tion factors used for models with local ISM dust content 
(D/G = 1). If this assumption is correct, a great increment 
of N gas abundance for Z/Zq > 1.5 should be observed, sug- 
gesting that it is not necessary to invoke an excess of N in 
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order to match the observed [N II] A6584 line intensity. It is 
important to stress that the depletion of the elements must 
be taken into consideration even if the nebula is without 
dust, such as in the radiation-pressure dominated models of 
I Groves et al.l (["2004) , according to which the dust is blown 
away by the wind. In our models, we assumed only two types 
of dust: graphite and silicates; their reference abundances 
are [C/H] = 1.22 • 10~ 4 when D/G = 1, in the graphite, 
and [O/H] = 1.94 • 10" 4 , [Mg/H] = 3.15 ■ 1(T 5 , [Si/H] = 
2.82 • 10" 5 , [Fe/H] = 2.70 ■ 10~ 5 for silicate, respectively, 
while the distribution of grain size is in the range [0.005, 
0.25] fan. 

The single cloud model is only the first step of our mod- 
cling effort, since in the observed Seyfert spectra two distinct 
components are usually detected, yielding high and low ion- 
ization lines. More accurate simulations can be produced in- 
troducing a double component model with two clouds, each 
accounting for either the high or the low ionization lines. We 
considered all the possible combinations between two clouds 
with fixed power law indices and metallicities. The spectral 
lines emitted by the first and second cloud are combined by 
a weighted mean, where the weight is given by the product 
of H/3 luminosity ratio and the ratio of the solid angles sub- 
tended by the clouds, as seen from the source of ionizing 
radiation {oji/uji, called geometrical factor and henceforth 
indicated with GF). The combined fluxes were derived from 
the geometry of the system. Given a line luminosity: 



: 4ivr 2 Sui/4 



(6) 



where r is the radius of the nebula, uj/4-k is the covering 
factor, S the emission line intensity (erg s _1 cm~ 2 ), from the 
relation between flux and luminosity, we get: 



(7) 



where F is the observed flux (ergs -1 cm -2 ) and d is our 
distance to the system. Considering two clouds, from the 
definition of ionization parameter in Eq. (3), we have that: 



Ui _ <3(H ) 47rr 2 n 2 c 
U 2 A-Krfnic Q(H ) 



so that: 



7j2 _ Urni 
r\ U 2 n 2 ' 

From Eq. (7) we also have the flux expressed as: 



Fi = 



4ird 2 



and 



■Fa _ i~2 S2W2 
which, recalling Eq. (9), finally gives: 

F2 _ U1TI1S2U2 

Fi L^r^SiWi 
In terms of H/3, this ratio can be written as: 
fa(H0) _ Qp £2 (H/3) 



(8) 



(9) 



(10) 



(11) 



(12) 



(13) 



Fi(H/3) — L\{B.py 
where L\ (H/3) /L* (H/3) is the H/3 luminosity ratio if wi = 0^2 



These quantities are calculated by CLOUDY for each model. 
If we introduce the emission line intensity 7a relative to H/3: 



F»(A) 

F tot {wy 



(14) 



where F to t (H/3) is the total observed flux of H/3, we obtain: 



hot{\) = 



+ Ji(A)F 1 (H/3)/F 2 (H/3) 
1 + Ji(H/3)/F 2 (H/3) 



and finally: 
/ to t(A) = 



MA) + h(X)GF ■ Li (H/3) /L 2 (H/3) 
l + GF-i*(H/3)/L*(H/3) 



(15) 



(16) 



If GF goes to zero, only the second cloud is visible; on the 
contrary, if GF goes to infinity, only the first cloud is visi- 
ble. Five values of GF have been used (0.25, 0.5, 1, 2 and 
4), so that each pair of models provides five new synthetic 
spectra. The total number of models obtained with the con- 
straints on D/G, a, Z/Zq, fixed for each pair of clouds, is 
approximately 9 ■ 10 6 . 



3.2.2 The P-index method 

The determination of chemical abundances in the H II re- 
gions of SFGs is complicated by the difficulty of getting re- 
liable measurements of the temperature of the H II clouds. 
The auroral lines, such as [O III] A4363, [S III] A6312 and 
[N II] A5755, used to evaluate the temperature, are usually 
weak, especially in low-excitation metal-rich H II regions. 
For this reason, the dir ect method i s almo st useless in the 
vast majority of cases. IPaeel et alJ l| 19791 ) suggested that 
some empirical calibrations between the nebular line fluxes 
could be used to derive directly the O abundance or the tem- 
perature of H II regions. Specifically, the parameter defined 



R23 = 



F([OII]A3727) + F([OIH]AA4959, 5007) 
FW) 



(17) 



was used to measure directly the abundance of O. This ra- 
tio was re-calibrated using photoionization models and con- 
firmed through other observations. Nonetheless, the com- 
parison between the results obtained from different samples 
of galaxies remains very difficult because of the high het- 
erogeneity of the model parameters and data. In all cases, 
the one-dim ensional calibrations are systematically wrong 
as shown bv lPilvugin fc Thuanl (|2005h . The R23 index alone 
does not remove all the degeneracies, so a more general two- 
dimensional parametric calibration, called the P-index, was 
proposed, introducing a new excitation parameter: 



P23 



J F([OIII]AA4959, 5007) 



(18) 



F([OII]A3727) + F([OIII]AA4959, 5007) ' 

IPilvugin fc Thuanl |2005l ) used their collection of high 
precision measurements to revise the P-index calibration, 
suggesting that the O abundance of H II regions in SFGs 
could be directly estimated, without degeneracies, from: 



12 + log(0/H) = 



(#23 + 726.1 + 842.2 ■ P + 337.5 ■ P 2 ) 
(85.96 + 82.76 ■ P + 43.98 • P 2 + 1.793 • #23) 



(19) 
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Kauffmann 

(M ) 

Figure 4. Comparison bet ween the average mass estimates de- 
rived by ISalim et all d200^ and the stellar masses computed by 
means of Eq. (20). The continuous line represents the perfect 
identity relationship, while the error bars describe the median un- 
ce rtainty estimates d erived from the spectral synthesis technique 
of ISalim et al,l l l2007f) and that descending from the photometric 
measurement errors. 



3.3 Stellar masses 

The correlations between the Super Massive Black Hole 
(SMBH) masses and the properties of the host stellar bulge, 
such as its luminosity, mass, and velocity dispersion, rep- 
resent the strongest argument in support for a common 
evolutionary track. Investigating these relationships in de- 
tail, however, requires large samples of high precision ob- 
servations, in order to properly calibrate empirical laws 
on the observed data. Combining the spectroscopic and 
photometric information , collected in the SDSS database, 
Kauffman n et all (|2003al ) developed an useful technique to 
estimate the stellar masses of AGN host galaxies. They ex- 
tensively exploited stellar population synthesis, in order to 
reproduce the observed colors and luminosities, thus deriv- 
ing the most appropriate mass-luminosity rati os for these 
object s. We adopted the updated results of ISalim et al.l 
l|2007l ) as the stellar masses of our S2G and SFG samples. 

The case of normal galaxies, however, is somewhat dif- 
ferent, since most of these objects are not included in any 
similar catalog, due to the absence of prominent emission 
lines in their nuclei. To estimate the stellar mass range 
covered by this fraction of our sample, we started from 
the results of our stellar population synthesis, which found 
no evidence of peculiarities in the stellar populations lo- 
cated in the circum-nuclear regions of NGs. From the stel- 
lar population synthesis of the spectra collected by the 
SDSS instrumentation, we derived the stellar masses in the 
circum-nu clear regions of these galaxi es and their character- 
istic M/L |Cid Fernandes et al.ll2005h . Then, assuming that 
M/L is constant throughout the galaxy, we computed the 
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Table 3. List of the emission lines calculated in the models 



Ion 


A (A) 


Ion 


A (A) 


ion 


A (A) 


r*~i tti 


3727 


[Ar 1 VJ 


4740 


[O 1J 


6363 


[O II] 


3729 


H/3 


4861 


[N II] 


6548 


[Ne III] 


3869 


[O III] 


4959 


Ha 


6563 


[Ne III] 


3968 


[O III] 


5007 


[N II] 


6584 


[S II] 


4064 


[N I] 


5200 


[S II] 


6716 


H7 


4340 


[Fe VII] 


5721 


[S II] 


6731 


[0 III] 


4363 


He I 


5876 


[Ar III] 


7135 


He II 


4686 


[Fe VII] 


6087 


[O II] 


7325 


[Ar IV] 


4711 


[OI] 


6300 







stellar masses from: 

M *°t = 10 °- 414 (''/.be,.-'4t) ) (20) 

where M t * ot and Alj iber are, respectively, the total mass in 
stars and the estimated stellar mass in the area covered by 
the fiber, while r' tot and r'j ibt , r represent the r' band Pet- 
rosian magnitude and the magnitude within the fiber. In 
principle, the uncertainty to be associated with this mass 
estimate results from two factors: at first, we have the photo- 
metric errors in the measurement of the fiber magnitude and 
in the assumption of the model representing the total galaxy 
flux; on the other hand there is the assumption of a constant 
representative M/L ratio across the whole galaxy. In order 
to test the reliability of our method, we compare the mass es- 
timates, deri ved with the photom etric uncertainty only, with 
the results of lSalim et alj l|2007l ). for the objects included in 
both samples (Fig. 4). Given the small uncertainty range of 
the photometric measurements, we expect that the the scat- 
ter in the correlation between the two methods is mostly a 
consequence of the assumption on M/L. The comparison 
of results gives a high degree of correlation (R — 0.978), 
which supports the use of the photometric measurements to 
provide an estimate of the total stellar masses. 



4 RESULTS AND DISCUSSION 

4.1 Comparison of models and observations 

To study the gas phase metallicities of active galaxies, we 
compared the observed emission lines with the synthetic 
spectra, produced by our models. In Table 3 we give a list 
of the lines that we included in the model calculations. In 
general we did not require that all the listed emission lines 
could actually be detected in the spectra, since some of them 
could be intrinsically weak or absent. In order to apply dif- 
ferent models to specific spectra, we only required that the 
number of detected emission lines was high enough to con- 
strain the parameter space of the models in use. We spent a 
considerable effort to reproduce as effectively as possible the 
observed spectra, with the ai m of determining the most re- 
alistic set of input parameters. lOliva. Marconi fe Moorwoodl 
I 19991 ) proposed a new method for deriving abundances 
in the NLR of active galaxies, consisting in a selection of 
fair models from a large set (27000 models, in their work), 
which is capable of fitting the observed line spectra, and 
then slightly modifying the chemical abundances in order to 
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determine the best estimates of the metallicity of the AGN. 
Our approach is slightly different. Instead of varying the 
input parameters, we produced a very large number of syn- 
thetic spectra, which we employed to determine the set of 
models, yielding the most accurate synthetic spectral fea- 
tures for each observed spectrum, with a x 2 test. In this 
case, the \ 2 ls given by: 

i 

where F° bs — F™° d is the difference between observed and 
modeled flux of the i-th line and at the associated error. 
Therefore, the most probable set of input parameters can 
be determined for each observed spectrum. The comparison 
has been carried out taking into account only the synthetic 
lines with a predicted intensity corresponding to a S/N > 3 
detection in the observed spectrum. Lines below the S/N 
observational threshold are excluded from the residual cal- 
culation. Another important concern, in the comparison of 
models with observed spectra, is the problem of intrinsic ex- 
tinction. While the observed emission lines are affected by 
a certain amount of reddening, the model spectra are ex- 
tinction free. Therefore, before performing the x 2 test, we 
artificially reddened our model spectra in order to achieve 
the observed F(Ha) / F(H/3) ratio. To do this, we combined 
F(Hq) ± AF(Hq) with F(H/3) ± AF(H(3) and we evaluated 
9 different Balmer line ratios on the observed spectra. As- 
suming a wavelength dependen t extinction law in the form 
relate d to the Av parameter (|Cardelli. Clayton fc Mathij 
1989), we derived all the extinction curves giving rise to a 
non-negative value of Ay ■ In this way, we produced up to 9 
different synthetic spectra for each model. 

The synthetic spectra are expressed in a relative flux 
scale, normalized to the predicted flux of H/3, F(Hf3), which 
acts as a scale factor. To compare these data with observa- 
tions, we converted the synthetic spectra in an absolute flux 
scale, by means of this factor. The accepted model has to 
satisfy a predetermined significance level, suggested by the 
analysis of the reduced \ 2 distribution, given by: 

Xred = ^ , (22) 

n, 7" Tip 

where n r is the number of measured lines and n p the number 
of parameters used in the models. In order to consider the 
models a good approximation of the observed spectra, the 
Xred distribution of the residuals left by the best fit mod- 
els over the sample of observations must tend to a Gaussian 
function, peaked at 1. A peak located at much smaller values 
would probably indicate an over-fitting, where non signifi- 
cant parameter sets fit very large uncertainty ranges. On 
the other hand, if the peak falls at larger values, the mod- 
els cannot reproduce the observed spectra, or the errors are 
too small. After several tests, we found that the best match 
between models and observations could be carried out as- 
suming the observed flux uncertainties in Eq. (21) at the 2a 
level. Since the model parameters affect each other in a com- 
plex interplay, we had to estimate the relative uncertainties 
by introducing a confidence range of 10 per cent around the 
best fit values. All the models falling in the resulting param- 
eter space region were taken into account and those meeting 
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Figure 5. The X^ed distribution for the D/G = 1 best fit mod- 
els carried out by single cloud (blue histogram), two clouds (red 
histogram), and Binette's models (black histogram) over the S2G 
galaxy sample. 



the condition: 

2 2^ Xmm /oo\ 

X - Xrmn < (23) 

(corresponding to a 2 5 per cent level of confidence, but see 
iMolla fc Hardy) l2002h were considered acceptable as well. 
The best fit parameters and their uncertainties are evaluated 
by computing the average and r.m.s. values from all the 
acceptable models. 

In our analysis, we consider 7 free parameters for the 
single cloud models and 11 for the two-cloud ones. As a con- 
sequence, more than 7 measured lines are required to fit the 
single cloud models, while at least 12 measured lines must 
be available for the composite models. About 50% of the 
spectra shows less than 12 measured lines, meaning that the 
two-cloud models can be applied only to 50% of the origi- 
nal galaxy sample. This fraction, however, includes spectra 
where S/N is high enough to detect also several weak lines 
above the 3-<j significance level. Single cloud models, on the 
other hand, have difficulties to reproduce with reasonable 
accuracy spectra with many observed emission lines. Since 
the simple rejection of spectra that we are not able to fit 
with some models would raise a strong bias on the inferred 
properties of our sample, we studied the different model ef- 
ficiencies in fitting the observed spectra. To achieve this, we 
split the S2G sample in two groups, namely spectra with 
n r > 12 and spectra with n r < 12. Since the D/G ratio is 
not considered as a model parameter, the models with dif- 
ferent D/G values are treated separately, so that we are left 
with four families of models (namely D/G — 0.25, 0.50. 0.75, 
and 1 times the Galactic reference). On a total of 2153, 1141 
spectra have n r < 12 and the remaining n r > 12. For this 
reason, the comparison between single cloud and two-cloud 
models is only possible for 1012 spectra. The total percent- 
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Figure 6. The cumulative \ 2 distribution for the D/G = 1 fam- 
ilies of single cloud (blue histogram), two cloud (red histogram), 
and Binctte's models (black histogram) applied to the shared 
spectra. 



age of spectra, which produced satisfactory fits, are ~ 65 
per cent and ~ 50 per cent with the single and two-cloud 
models respectively. The percentage of spectra successfully 
fitted with n r < 12 is 80 per cent, while ~ 45 per cent is 
the fraction of spectra fitted with n r > 12. The Xred distri- 
butions for the D/G = 1 model families are shown in Fig 5. 
The distributions of the single cloud and Binette's models 
for high D/G peak at Xred ~ 0.6, while the two-cloud mod- 
els show tails of high residuals, although their peak is closer 
to the desired value of Xred = !• 

A more qualitative difference between the model fam- 
ilies becomes evident by comparing the distribution of ab- 
solute x 2 illustrated in Fig. 6, since this statistics repre- 
sents the goodness of the overall fitting procedure (the in- 
put model parameters, the choice of the significance level 
and the errors) and it is possible to discriminate between 
two models able to fit the same spectrum. For each assumed 
D/G value, there are respectively 289, 356, 410, and 396 
spectra, which were fitted by single cloud, two-cloud and 
Binette's models, but with dramatically different x 2 distri- 
butions. The average values of Xred f° r each group, with 
the corresponding standard deviations, are listed in Table 4, 
but the differences can be better appreciated looking at the 
cumulative distributions of x 2 plotted in Fig. 6: when ap- 
plied to the same spectra, two-cloud models produce lower 
residuals than single cloud or Binette's models. This is rea- 
sonable because spectra with a larger number of measured 
lines could be more effectively fitted with an increased pa- 
rameter resolution, particularly for the power law index and 
the ionization structure, at the cost of a higher computa- 
tional load. Anyway, these considerations should be cau- 
tiously handled, since the percentage of successful fits are 
calculated averaging over few spectra, so that a single mea- 
surement could drastically modify the result. The most af- 



Table 4. Comparison of the average residuals left by single and 
two-cloud models on the shared spectra. 



D/G 


-^Vspectra 


(X re( ^)single cloud 


{X re£ ^) two — cloud 


0.25 


289 


0.929 ± 0.373 


1.207 ± 0.590 


0.50 


356 


0.914 ± 0.398 


1.182 ± 0.636 


0.75 


410 


0.922 ± 0.418 


1.140 ± 0.607 


1.00 


396 


0.862 ± 0.407 


1.112 ± 0.578 



fected lines are those from [Ar IV] AA4711, 4740, [N I] A5200, 
[Fe VII] AA5721,6087, He I A5876, and [Ar III] A7135. 



4.2 The properties of gas 

In Table 5 we summarize all the mean parameter values giv- 
ing the best fit solutions, with their la uncertainty ranges, 
arranged according to the different model families. We ob- 
serve that the power-law index a of the ionizing spectrum 
decreases to -1.8, as soon as we increase the D/G ratio, 
probably because of the decreasing gas metallicity content 
(a large fraction of metals is depleted into grains). Indeed, a 
low metallicity of the gas implies a much slower cooling and 
a weaker ionizing spectrum can produce the same degree of 
ionization. The total metallicity is well defined in all models 
with negligible scatter: Z — 1.5 — 2.5 Zq is the range of fit- 
ted values with the single cloud models, while the two-cloud 
models point towards estimates of Z = 1.0 — 2.0 Zq. The 
high metallicity of the gas indicates that the AGN fuel is an 
evolved material. This, in turn, supports the idea that the 
gas has a local origin in the vast majority of the observed 
S2Gs. The densities are in agreement with the expected val- 
ues, since single cloud and two-cloud models produce sim- 
ilar density estimates, both consistent with the indications 
obtained from the [S II] AA6716, 6731 intensity ratios. The 
column densities are very similar in the single cloud mod- 
els and in the first cloud of the two-cloud models, while the 
second cloud of the two-cloud models has a lower column 
density and a lower ionization parameter, as well, allowing 
us to treat such cloud as an ionization bounded case. The 
ionization parameter in the single cloud models is almost 
constantly fixed at log((7) = —3 for different values of D/G. 
In two-cloud models the high ionization cloud shows a little 
variation of the ionization parameter, depending on D/G, 
with a mean value of log(?7) = —2.6, while very steady val- 
ues of log(C/) = —3.3 are found for the low ionization cloud. 
In general, the amount of ionization in the single cloud mod- 
els is intermediate between the values found in the first and 
second cloud of the two-cloud models. 

On the other hand, as we show in Fig. 7, a substan- 
tially different situation applies to the chemical properties 
of the ISM in star forming galaxies. According to the metal- 
licity indicators discussed in §3.2.2, these objects are rather 
characterized by ISM with a sub-solar chemical composi- 
tion, suggesting that their on-going star formation is a first 
generation process. This result is also confirmed by the ev- 
idence that stellar masses in galaxies hosting nuclear star 
formation are gene rally lower than in S2Gs. in agreement 
with the findings of iKauffmann et all l|2003bl ). 
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Table 5. Average output parameters and standard deviations for the four D/G families of single cloud, two-cloud and 
Binette's photoionization models. 



Parameter 


SCo.25 


SC0.50 


SCo.75 


SCi.oo 


2Co.25 


2Co.50 


2Co.75 


2Ci.oo 


Bincttc 


a 


-1.5 ± 0.3 


-1.6 ± 0.3 


-1.7 ± 0.3 


-1.8 ± 0.3 


-1.3 ± 0.3 


-1.4 ± 0.3 


-1.6 ± 0.3 


-1.8 ± 0.3 


-2.0 ± 0.2 


z/z e 


2.0 ± 0.6 


1.9 ± 0.6 


1.9 ± 0.5 


1.8 ± 0.5 


1.6 ± 0.5 


1.7 ± 0.4 


1.7 ± 0.4 


1.7 ± 0.4 


1.5 ± 0.4 


logn e ,i 


2.3 ± 0.7 


2.3 ± 0.7 


2.2 ± 0.7 


2.1 ± 0.7 


2.2 ± 0.8 


2.2 ± 0.8 


2.1 ± 0.7 


2.2 ± 0.8 


n.a. 


log7V c ,i(H II) 


20.5 ± 0.7 


20.3 ± 0.6 


20.0 ± 0.5 


19.9 ± 0.4 


20.2 ± 0.4 


20.2 ± 0.4 


20.2 ± 0.4 


20.2 ± 0.4 


20.0 ± 0.4 


logf/i 


-3.1 ± 0.2 


-3.1 ± 0.2 


-3.0 ± 0.2 


-3.0 ± 0.3 


-2.7 ± 0.4 


-2.6 ± 0.4 


-2.6 ± 0.4 


-2.5 ± 0.4 


-1.6 ± 0.2 


logn e ,2 


n.a. 


n.a. 


n.a. 


n.a. 


3.0 ± 0.9 


2.8 ± 0.9 


2.7 ± 0.9 


2.5 ± 0.8 


n.a. 


log7V c , 2 (H II) 


n.a. 


n.a. 


n.a. 


n.a. 


19.7 ± 0.4 


19.7 ± 0.3 


19.7 ± 0.3 


19.7 ± 0.3 


n.a. 


logU 2 


n.a. 


n.a. 


n.a. 


n.a. 


-3.3 ± 0.2 


-3.3 ± 0.2 


-3.3 ± 0.2 


-3.3 ± 0.3 


-3.5 ± 0.3 


GF 


n.a. 


n.a. 


n.a. 


n.a. 


n.a. 


n.a. 


n.a. 


n.a. 


4.1 ± 2.1 


A v 


1.4 ± 0.6 


1.5 ± 0.6 


1.5 ± 0.6 


1.6 ± 0.6 


1.2 ± 0.4 


1.2 ± 0.4 


1.3 ± 0.4 


1.4 ± 0.4 


1.6 ± 0.6 
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Figure 7. Relationship between the estimated stellar masses and Figure 8. Metallicity of the circum- nuclear stellar population, 

gas phase metallicitics for SFGs (blue points) and S2Gs (red plotted as a function of the population age of SFGs (blue points), 



points). 



S2Gs (red points), and NGs (black points). 



4.3 The properties of stars 

The processes of star formation and stellar feedback origi- 
nate an interplay, which connects the chemical properties of 
stars with those of the ISM, wherein they reside. The inter- 
stellar gas is chemically enriched through stellar activity, 
while the metallicity of stars reflects the chemical composi- 
tion of the environment, where they are formed. The spec- 
tral synthesis technique, which we exploited to account for 
the stellar continua and absorption features, reproduce the 
circum-nuclear stellar populations of o ur galaxy sample, in 
terms of a mixture of SSPs. Following ICid Fernandes et al.l 
(2005), we estimate the average stellar population metallic- 
ities through: 

92 

(Z*) = ^ Mi Z i , (24) 

i=l 



where fii is the mass fraction ascribed to the i-th SSP com- 
ponent of the spectrum. Considerations on the application 
of such paradigm to large samples led to evaluate a typical 
uncertainty of AlogZ* ~ O.ldex. 

Despite the different chemical compositions of the gas 
seen in Fig. 7, the average metallicity detected in the stel- 
lar populations in SFGs and S2Gs is mostly similar, but 
with a slight over-abundance of heavy elements in the sec- 
ond group. As it is shown in Fig. 8, indeed, the metallicity 
of stars in the SFG sample stays low, with no appreciable 
dependence on the age of the populations, as we would ex- 
pect for a generation of stars originated by an unpolluted 
medium. S2Gs, on the contrary, exhibit also a spread of 
high metallicity population components, which represents 
a hint of further star formation in a chemically enriched 
medium, bringing such objects on the same range of metal- 
licities, commonly observed in the circum-nuclear regions of 
NGs. This result indicates that SFGs and S2Gs host stel- 
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Figure 9. Relationship between the age of the nuclear stellar 
population and the total stellar mass found in the SFG, S2G, 
and NG samples. Symbols are the same as in Fig. 8. 



Figure 10. D n (4000A) spectral index plotted against the total 
estimated mass of SFGs, S2Gs, and NGs. Symbols are the same 
as in Fig. 8. 



lar populations mainly formed in a pristine inter-stellar gas, 
but, as it is also suggested by the average spectral proper- 
ties illustrated in Fig. 3, the process occurred earlier in S2Gs, 
where the O/B spectral class stars have already evolved off 
the main sequence. If this scenario is actually correct and we 
do observe stellar populations formed by pristine gas both 
in SFGs and S2Gs, then it is likely that the circum-nuclear 
regions of active galaxies experienced past star formation, 
on a time scale which should account for the evolution of 
massive and short-lived O/B stars and for the enrichment 
of the ISM, but not for that of lower mass A stars, as well 
as for the formation of new stars in the enriched medium. 

A look at the relationship between the stellar popula- 
tion ages in the circum-nuclear regions of galaxies and their 
total stellar masses, illustrated in Fig. 9, also reveals some 
further interesting effects. SFGs and S2Gs, indeed, lie on 
what appears to be a well defined sequence of increasing 
nuclear population age for larger total mass. The sequence 
has a high mass break and the range of masses, where AGN 
activity is observed, is also populated by NGs. The total 
stellar mass, therefore, increases regularly with age, through 
the star forming and nuclear activity stages. Later on, the 
stellar mass assembly is stopped and galaxies evolve as quies- 
cent objects, dominated by old populations. Before drawing 
any conclusions, however, we have to consider the possible 
biases, introduced by our sample selection criteria. In par- 
ticular, the choice of a S/N limiting threshold reflects in 
the selection of more luminous objects. SFGs are not par- 
ticularly prone to this problem, because the high luminosity 
of their young stellar populations enables us to detect ob- 
jects with fairly small masses, but for S2Gs and NGs this 
could be a real issue. If we look at the mass distribution 
of the SDSS line emitting sources, which, in the considered 
redshift range, fall in the AGN region of the O123 diagram, 
however, it turns out that a S/N > 10 threshold is not sys- 



tematically excluding low mass objects from the sample, so 
that the observed mass range for S2Gs is a real feature in 
the considered volume. The NG sample, on the other hand, 
experiences a quite different situation. In this case, we are 
not able to include very low mass quiescent objects, but the 
issue is not particularly critical for our interests, because the 
artificial cut-off falls below the transition region between star 
forming and active galaxies. Moreover, as we noted in the 
analysis of the average spectral properties, in order for the 
difference between the spectra of S2Gs and NGs to have an 
appreciable physical meaning, we need that the NG sam- 
ple covers a similar mass range, with respect to the S2Gs. 
Therefore, the adopted limiting constraints do not affect the 
S2G distribution, though we may expect that the NG popu- 
lation extends well below the detected lower end of the S2G 
mass range. Further investigation on the properties of low 
mass NGs would represent a useful extension of our work, 
but the task would require extensive spectral modeling of 
large numbers of spectra, which is beyond the scope of our 
current analysis. 

The differences that we appreciated on the basis of our 
averaged analysis can be better understood by looking at 
the diagram plotted in Fig. 10, where we show the dis- 
tribution of the D n (400 A) spectral index (defined as in 
iKauffmann et ail 1200 3a) as a function of the total stellar 
mass. The region covered by SFGs and S2Gs represents a 
well defined sequence of regular growth in the index, with 
increasing host galaxy stellar mass. Quiescent galaxies, in- 
stead, are well separated from this sequence. The situation 
is fully consistent with a smooth transition of gradually in- 
creasing stellar ages in a generation of stars formed in pris- 
tine gas, as soon as we move from low mass SFGs to high 
mass SFGs, all the way across S2Gs. The transition region 
can be identified within a range of total stellar masses of 
9.8 < log(M*/M0) < 10.3. The overlapping mass range 
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may represent a region where star forming activity is most 
likely going to start an interplay with the AGN. Depending 
on the interaction between these processes, the galaxy may 
either stop its stellar mass assembly, eventually evolving as a 
low mass quiescent system, or form new generations of stars, 
like the high metallicity populations detected in some S2G 
spectra of our sample, to become a high mass system. The 
difference between the two evolutionary possibilities proba- 
bly resides in the power of the nuclear source and the way 
it affects ISM in the surrounding space. 



5 CONCLUSIONS 

As we subdivide our observational sample in objects with 
different types of nuclear activity, taking into account 
Seyfert 2, star forming nuclei and normal galaxies, we ob- 
serve that the circum-nuclear environments are character- 
ized by appreciably different properties. Objects featuring 
ionized gas in their nuclear spectra show distinct ISM chem- 
ical compositions, characteristic of pristine gas in SFGs and 
of evolved medium in S2Gs. The stellar populations, how- 
ever, do not reflect such a dichotomy, being almost every- 
where consistent with star formation powered by pristine 
gas. Another striking property is the appreciably tight rela- 
tionship, which holds between the age of stars and the global 
stellar mass of galaxies. A relationship between the age of 
stars and the stellar mass observed in galaxies is expected, if 
more massive systems are also forming more quickly. Since 
the circum-nuclear regions of active galaxies appear to be 
more likely explained as the result of a single star forma- 
tion event, rather than as the combination of several loops, 
our analysis suggests that nuclear activity takes place after 
the bulk of stellar mass assembly has occurred. Because of 
its role in controlling the star formation time-scale, the ob- 
served stellar mass is likely a key parameter in determining 
whether an AGN can be triggered and if star formation can 
affect the circum-nuclear regions further on. 

The average trends that we identified in our sample ac- 
tually mark an important distinction in the environment of 
individual galaxies, whose nuclei are affected by different 
degrees of activity. The properties observed within a wide 
sample of galaxies, including normal objects, as well as star 
forming and active ones, are consistent with the expectation 
that star formation and nuclear activity should be actually 
related in the evolution of galaxies. The connection can be 
largely supported on the basis of sound theoretical grounds, 
but it is still missing a conclusive observational confirma- 
tion. In our analysis we identified several properties of the 
selected galaxies, which may lead to a more precise investi- 
gation of the problem and, ultimately, to a clearer view of 
the general picture. 

The main results, which we are able to point out on the 
basis of our data, appear as follows: 

(i) there is a range of galaxy masses (9.8 < 
log(M»/M©) < 10.3) where SFGs and AGNs are character- 
ized by circum-nuclear stellar populations of approximately 
the same age and similar chemical abundances in the stellar 
component; 

(ii) objects in this range are the oldest SFGs and the 
youngest AGNs (in terms of stellar population); 



(iii) although the metallicity of stars in SFGs within this 
mass range are indistinguishable from those of lower mass 
objects and from most of the S2Gs, there are some S2Gs 
showing evidence for further evolution. 

The metallicity of the stars seems to be independent 
from the mass of the whole galaxy and from the age of 
the circum-nuclear stellar population, as it is the case for 
the gaseous component. This suggests that we are observing 
different stages of the evolution of the same star formation 
process. Therefore, the circum-nuclear regions of SFGs with 
mass 9.8 < log(M*/MQ) < 10.3, which are clearly in a more 
advanced evolutionary stage than lower mass SFGs, seem 
to be the most promising candidate precursors of AGNs, if 
there is any causal link between the two phenomena. These 
results are in good agreement with the conclusions reported 
in the high resolution studies of lDavies et al.l l|2007h . who, in- 
vestigating the nuclear environment on a much smaller scale, 
in a few particular cases, were able to detect co-existing star 
formation and nuclear activity and to estimate the typical 
time-scales connecting the two processes. Further evidence 
for a connection between star formation and nuclear activ- 
ity from a much more similar dataset was also achieved by 
IShavinski et al.l 120071 ) . in their analysis of the role of AGN 
feedback in the evolution of galaxies. Looking at a sample 
of early type objects in the redshift range 0.05 < z < 0.1, 
they concluded that AGN feedback is a key process in the 
evolution of a galaxy through its stages of star formation, 
nuclear activity and quiescent existence. They also conclude 
that the AGN feedback process is likely to smoothly affect 
the star formation rate of a galaxy, during time, and that 
more massive objects experienced this transition earlier. 
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